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Good geometrical calibration is essential in the use of high resolution detectors. The individual 
sensors in the detector have to be calibrated with an accuracy better than the intrinsic resolution, 
which typically is of the order of 10 /^m. We present an effective method to perform fine calibration 
of sensor positions in a detector assembly consisting of a large number of pixel and strip sensors. 
Up to six geometric parameters, three for location and three for orientation, can be computed for 
each sensor on a basis of particle trajectories traversing the detector system. The performance of 
the method is demonstrated with both simulated tracks and tracks reconstructed from experimental 
data. We also present a brief review of other alignment methods reported in the literature. 



I. INTRODUCTION 

For full exploitation of high resolution position sen- 
sitive detectors, it is crucial to determine the detec- 
tor location and orientation to a precision better than 
their intrinsic resolution. It is a very demanding task 
to assemble a large number of detector units in a large 
and complex detector system to this high precision. 
Also, after assembly, the position determination of the 
modules by optical survey has its limitations because 
of detectors obscuring each other. Therefore the fi- 
nal tuning of detector and sensor positions is made by 
using reconstructed tracks. 

In this paper we present an effective method by 
which individual sensors in a detector setup can be 
aligned to a high precision with respect to each other. 
The basic idea is illustrated in Figure ^ Using a large 
number of tracks, an optimum of each sensor position 
and orientation is determined such that the track fit 
residuals are minimized. 
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FIG. 1: Schematic illustration of the method to correct 
the sensor position: using a large number of tracks i,i = 
1,...,N and measured hits on a detector, the sensor is 
moved such as to minimize the residuals. 



Sections and IVII we demonstrate the performance 
of the method applied to a test beam setup and to 
a simulated pixel vertex detector, respectively. The 
CMS % Pixel detector is is used as a model in the 
simulation. 



II. BRIEF REVIEW OF ALIGNMENT 
METHODS 

Most HEP experiments equipped with precise 
tracking detectors have to deal with misalignment is- 
sues, and several different approaches for alignment 
by tracks have been used and reported. Most meth- 
ods are iterative with 5-6 parameters solved at a time. 

Several papers concerning different aspects of align- 
ment in the DELPHI experiment can be found in the 
literature. For instance, — ^ M^M" ^^'^ cosmic 
rays are used for the global alignment between sub- 
detectors VD, OD and TPC |||. The most detailed 
DELPHI alignment paper deals with the alignment of 
the Microvertex detector Q . 

In the ALEPH experiment, alignment was carried 
out wafer by wafer, and with 20 iterations and 20000 
Z° — > qq and 4000 Z° — > events an accuracy of 

a few fj,m can be achieved [3- 

A different, computationally challenging approach 
is chosen in the SLD experiment, where the algorithm 
requires simultaneous solution of 576 parameters lead- 
ing to a 576 by 576 matrix inversion In the SLD 
vertex detector, a recently developed matrix singular 
value decomposition technique is also used for internal 
ahgnment @- 



III. COORDINATE SYSTEMS AND 
TRANSFORMATIONS 



The outline of this paper is as follows: In Section 
^we briefly review published alignment methods. In 
Section [nil we introduce the basic notations and coor- 
dinate systems involved in our method. In Section ITVI 
we present the detailed formulation of the method. In 



Our method is applicable to detector setups which 
consist of planar sensors like silicon pixel or strip 
detectors. For track reconstruction one convention- 
ally uses the local (sensor) coordinate system and the 
global detector system. The local system {u,v,w) is 
defined with respect to a detector module (sensor) as 
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follows: The origin is at the center of the sensor, the 
u)-axis is normal to the sensor, the w-axis is along 
the precise coordinate and the u-axis along the coarse 
coordinate. The global coordinates are denoted as 
{x,y,z)- 

The transformation from the global to the local sys- 
tem goes as: 



q = R(r - ro) 



(1) 



where r = {x,y,z), q = [u, v, w), R is a rotation and 
To = (a^OjJ/Oi^o) is the position of the detector center 
in global coordinates. 

In the very beginning of the experiment the rota- 
tion R and the position rg are determined by detec- 
tor assembly and survey information. In the course of 
experiment this information will be corrected by an 
incremental rotation AR and translation Ar so that 
the new rotation and translation become: 



R 



ARR 

ro + Ar. 



The correction matrix AR is expressed as: 



AR 



(2) 
(3) 



(4) 



where Rq , R^ and R^ are small rotations by 
Aa, A/3, A7 around the w-axis, the (new) w-axis and 
the (new) ly-axis, respectively. The position correc- 
tion Ar transforms to the local system as: 



Aq = ARRAr 



(5) 



with Aq — {Au, Av, Aw). Using (|1I5|I we find the 
corrected transformation from global to local system 
as: 



q= = ARR(r ro) Aq. 



(6) 



where the superscript c stands for 'corrected'. The 
task of the alignment procedure by tracks is to de- 
termine the corrective rotation AR and translation 
Ar or Aq as precisely as possible for each individual 
detector element. 



The equation of a straight line in global coordinates, 
approximating the trajectory in a vicinity of the de- 
tector, can be written as: 



rs{h) =rx + hs 



(7) 



where is the trajectory impact point on the detec- 
tor in question, s is a unit vector parallel to the line 
and /i is a parameter. Equation 10 is for uncorrected 
detector positions. 

Using Eq. ^ the corrected straight line equation 
in the local system reads: 



qs{h) = Rc(ra; + h s - Tq) - Aq 



(8) 



where Rc = ARR. A point qs — qs(/ix) which lies 
in the detector plane must fulfill the condition q^ • 
w = 0, where w = (0, 0, 1) is normal to the detector. 
From this condition we can solve the parameter 
which gives the corrected impact or x-ing point on the 
detector: 



[Aq - Rc(ra; - ro)] ■ w 

RrS • W 



(9) 



The corrected impact point coordinates q^ in the local 
system are then: 

c T} f N , [Aq-Rc(r:r -ro)] • W ^ 

q^ = Rc{rx - ro) H ^ — RcS - Aq. 

RcS • w 

(10) 

Since the uncorrected impact point is q^; = R(rj^— ro), 
Eq. 110() can be written as: 

q-^ARq^+^^q-^^q-^-^ARt Aq. (11) 
ARt • w 

where t — Rs is the uncorrected trajectory direction 
in the detectors local frame of reference. Eq. 
evaluates to: 



q;^ = ARq^, + (Au; - [ARq^jg) 



ARt 

[ARt]s 



-Aq. (12) 



This expression provides us with a 'handle' by which 
the unknowns Aq and AR can be estimated by mini- 
mizing a respective function using a large number 
of tracks. 



IV. DESCRIPTION OF THE ALIGNMENT 
ALGORITHM 

A. Basic Formulation 

Since the alignment corrections are small, the fitted 
trajectories can be approximated with a straight line 
in a vicinity of the detector plane. The size of this 
small region is determined by the alignment uncer- 
tainty which is expected to be at most a few hundred 
microns so that the straight line approximation is per- 
fectly valid. 



B. General Solution 

We denote a measured point in local coordinates 
as q™ = {um,Vm,0)- The corresponding trajectory 
impact point is q!^ = {ux,Vx,0). For simplicity we 
omit the superscripts c in the coordinates and Vx- 
In stereo and pixel detectors we have two measure- 
ments. Urn and v„i, and in non-stereo strip detectors 
only one, In the latter case the coarse coordinate 
Vm is redundant. The residual is either a 2-vector: 
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or a scalar e — Eu — ~ Um- In the following we 
treat the more general 2-vector case. The scalar case 
is a straightforward specification of the 2-vector for- 
malism. 

The function to be minimized for a given detec- 
tor is: 



(14) 



where the sum is taken over the tracks j. Vj is 
the covariance matrix of the measurements {urn,Vm) 
associated with the track j. The alignment correc- 
tion coefficients, i.e. the three position parameters 
(Au, Av, Aw) and the three orientation parameters 
(Aa,A/3, A7) are found iteratively by a general 
minimization procedure. At each step of the itera- 
tion one uses the so far best estimate of the alignment 
parameters in the track fit. 

Let us denote these parameters as p = 
(Au, Aw, Aui, Aa, A/3, A7). Then, according to the 
general solution, the iterative correction to p has 
the following expression: 

-1 



Sp 



E 



E 



(15) 



where 3j is a Jacobian matrix of ej (p) : 

J, =Vp£,(p). (16) 

An adequate starting point for the iteration is a null 
correction vector p=0. 

In the general case of two measurements {um, Vm), 
Jj is a 6 X 2 matrix. In case of scalar e, for single 
sided strip detectors, Jj is a vector of 5 elements, be- 
cause Av is redundant and cannot be fitted. It will 
also be foreseen that only a sub-set of the 6 alignment 
parameters would be fitted and the others kept fixed. 
In this case the dimension of the Jacobian matrix re- 
duces accordingly. 

The derivatives of the Jacobian matrix can be com- 
puted to a good precision in the small correction angle 
approximation (see below) . The elements of the ma- 
trix J for a given track are then: 

f -I \ 

-1 
tan ip tan 

Vxtamp Wa;tani9 

Uxtanip Uxtan'd 

\ -Ux J 



(17) 



The quantities tan -0 and tan ■& are defined in the next 
section. 



C. Linearized Solution with the Tilt Formalism 

We call "tilts" the angle corrections x which are 
small enough to justify the approximations cos a: ~ 1 



and sin a; ~ x. In this approximation the correction 
matrix AR reads: 



1 A7 A/3 \ 
AR= I -A7 1 Aa 
-A/3 -Aa 1 J 



(18) 



Using Eq. lfTF|l we linearize Eq. ifT^ and get the fol- 
lowing expressions for the corrections of the impact 
point coordinates as a function of the alignment cor- 
rection parameters: 

Aux = — Au + 6 ta.n ijj + A'y Vx (19) 
Avx = -Av + Stan-d ~ A^Ux (20) 

where S = Aw -t- A/3 Ux + Aa Vx ■ The quantity ip is 
the angle between the track and the uw— plane and 
■d is the angle between the track and the uw— plane: 
tantp = h/h, tani? = ^2/^3- 

With this approximation the residuals (|13|l depend 
linearly on all 6 parameters. Hence the minimiza- 
tion problem is linear and can be solved by standard 
techniques without iteration. 



V. ALIGNMENT OF A TEST BEAM SETUP 

From Eqs. (|19|l and H2U|I we can estimate the contri- 
butions of various misalignments to the hit measure- 
ment errors. For example the contribution of a mis- 
alignment Aa around the u-axis to the w-coordinate 
is: 



Aw ~ V Aatan-d. 



(21) 



The error is small near normal incident angles, but 
grows rapidly as a function of i?. At i? = 45° and near 
the edge of the sensor {v = 3 cm) the error goes as 
30000 fim Aa so that for only 1 mrad error in Aa the 
systematic error in the w-coordinate is 30 /im. 




FIG. 2: 
beam. 



Helsinki Si Beam Telescope in the CERN H2 



The silicon detector team of Helsinki Institute of 
Physics made a precision survey of detector resolution 
as a function of the angle of incidence of the tracks . 



TULT008 



4 



CHEP03, La Jolla California, March 24-28, 2003 



The study was made in the CERN H2 particle beam 
with a setup described in Figure [3 One of the sihcon 
strip detectors was fixed on a rotative support which 
allowed the tracks to enter between and 30 degrees 
of incident angle. The angular dispersion of the beam 
was about lOmrad and the hits covered the full area 
of the test detector. 

In order to obtain reliable results it was extremely 
important to calibrate the tilt angle to a very high 
precision. Our algorithm was used in the alignment 
calibration. In Table^we show the result of the align- 
ment demonstrating the precision obtained by about 
3000 beam tracks. 

TABLE L Alignment parameters obtained by the algo- 
rithm 



Parameter 


At degrees 


At 30 degrees 


A?i(/im) 
Aw(pm) 
AQ(mrad) 
A/3(mrad) 
A7(mrad) 


186.0±0.1 
200±20 
5.6±0.7 
5.8±0.9 
-14.12±0.01 


-264.7±0.1 
-131±6 
12.9±0.9 
32.59±0.04 
-15.86±0.01 



With the precise alignment we have been able to 
determine the optimal track incident angle which min- 
imizes the detector resolution 8l. 



4. Reconstruction of the particle trajectories in 
the nominal (ideal) detector thus using slightly 
wrong hit positions. This simulates the realistic 
situation in which the detector alignment is not 
yet performed. 

For the simulated detector type we choose a vertex 
detector which is a simplification of the CMS Pixel 
barrel detector P, |lO| with two layers. The setup is 
illustrated in Figure There are 144 sensors in layer 
1 and 240 sensors in layer 2. The distance of the layer 
1 from the beam line is about 4 cm and the layer 2 
about 8 cm. 




FIG. 3: Illustration of the simulated vertex detector in the 
alignment study. 



VI. MONTE CARLO SIMULATION 

A. Simulated Detector 

A Monte Carlo simulation code was written to test 
the alignment algorithm. High momentum tracks 
were simulated and driven through a set of detector 
planes. The simulated hits were fluctuated randomly 
to simulate measurement errors. Gaussian multiple 
scattering was added quadratically using the Highland 
approximation. The algorithm involves misalign- 
ment of a detector setup in order to simulate a realistic 
detector. 

The experimenters' imperfect knowledge of the true 
position of the detector planes is simulated by recon- 
structing the trajectories in the ideal (not misaligned) 
detector. This means that in the transformation from 
local to global coordinate system one uses the ideal 
positions of the detector planes. The full algorithm in 
brief is as follows: 

1. Creation of an ideal detector setup with no mis- 
alignments 

2. Creation of a misaligned, realistic detector 

3. Generation of the particles and hits in the mis- 
aligned detector simulating the real detector 



In the simulation we used the following conditions: 

1. Misalignment of chosen sensors: The shifts 
Au, Au, Aw were chosen at random, each in the 
range ±100 /xm and the tilts Aa, A/3, A7 were 
chosen at random each in the range ±20mrad. 

2. Beam and vertex constraints: The vertex posi- 
tions were Gaussian fluctuated around the cen- 
ter of the beam diamond with Ux = ^ 20 /im 
and Uz = 7 cm and the tracks were fitted with 
the constraint to start from one point, i.e. from 
the primary vertex. 

In the following we consider two different cases of mis- 
ahgned detectors: 

I. All sensors in layer 2 fixed, all sensors in layer 1 
misaligned. 

II. Only one sensor in layer 2 fixed, all remaining 
383 sensors misaligned. 

In case I the total number of fitted parameters is 
6 X 144 = 864 and we used about 2 x 10^ tracks. 
The case I appears to be an 'easy' one with which 
the algorithm copes very well, as we see below. The 
second case we call 'extreme' since the alignment is 
based on one reference sensor which covers only about 
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Iteration of 6 alignment parameters for a Pixel sensor 
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Iteration of 6 alignment parameters for a Pixel sensor 
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FIG. 4: The six plots on the left show the rate of convergence in case I for the alignment parameters of one sensor 
(circular spots). The solid line shows the true parameter value. The six plots on the right are for the case II. 



0.26% of the detector setup area. The total number 
of fitted parameters in this case was 6 x 383 = 2298. 
In the following sectios we show perfomance results of 
the algorithm in these two cases. 



B. Convergence of the Algorithm 

The convergence rate of the alignment procedure as 
a function of the iteration cycle is shown in Figure 0| 
It appears that the convergence is fast in the 'easy' 
case (the 6 plots on the left) where more than 60% 
of the sensors provide the reference. The convergence 
takes place after a couple of iterations. 

In the case where only one sensor is taken as a ref- 
erence (plots on the right of the figure), the situation 
is different. It appears that the number of iterations 
needed varies between 20 and 100 from parameter to 
parameter. It is also seen that the converged parame- 
ter values are somewhat off from the true values, but 
the precision is reasonable. 



C. Comparison of Fitted and True Parameters 

The precision of the fitted parameters in compari- 
son with the true values is shown in Figure [S] on the 
left for the case I. The correlations are very strong. 



The typical deviation of the fitted parameters from 
the true value is less than 1 /im for the offsets and 
a fraction of a milliradian for the tilts. The preci- 
sion appears to be better than actually needed in this 
case, indicating that a smaller statistics would give a 
satisfactory result. 

In case II (the plots on the right of the figure) a good 
correlation is observed, but the precision is somewhat 
more modest. For example the error in Aw (the shift 
normal to the sensor plane) is still in most cases below 
10 /zm. 



VII. CONCLUSIONS 

We have developed a sensor alignment algorithm 
which is mathematically and computationally simple. 
It is based on repeated track fitting and residuals op- 
timization by minimization. The computation is 
simple, because the solution involves matrices whose 
dimension is at most 6x6. The method is capable of 
solving simultaneously all six alignment parameters 
per sensor for a detector setup with a large number of 
sensors. 

We have successfully applied the method in a 
precision survey of silicon strip detector resolution 
as a function of the tracks incident angle. Further- 
more, we have demonstrated the performance of 
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Fitted versus true parameters - all misaligned sensors 
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Fitted versus true parameters - all misaligned sensors 
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FIG. 5: The scatter plot of fitted versus true alignment parameters. There are 6 x 144 entries in the plots on the left 
(case I) and 6 x 383 entries in the plots on the right (case II). 



the algorithm in case of a simulated two-layer pixel 
barrel vertex detector. The method performs very 
well in the case where the outer layer is taken as 
a reference and all inner sensors are to be aligned. 
The algorithm performs reasonably well also in the 
extreme case where only one sensor, representing 
some 0.26% of the total area, is taken as a reference 
for the alignment. 
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